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O ! (Dated: March 15, 2012) 

■ We discuss the subleading contact interactions, or counterterms, of the triplet channels 
of nucleon-nucleon scattering in the framework of chiral effective field theory, with S and 
P waves as the examples. The triplet channels are special in that they allow the singu- 
lar attraction of one-pion exchange to modify Weinberg's original power counting (WPC) 
scheme. With renormalization group invariance as the constraint, our power counting for 
the triplet channels can be summarized as a modified version of naive dimensional analysis 

^ in which, when compared with WPC, all of the counterterms in a given partial wave (leading 

or subleading) are enhanced by the same amount. More specifically, this means that WPC 

■ needs no modification in 3 Si — 3 Di and 3 Pi, whereas a two-order enhancement is necessary 

^ ■ in both 3 P and 3 P 2 - 3 F 2 . 

ON 
ON 
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I. INTRODUCTION 

In applying the idea of chiral effective field theory (EFT) to nuclear physics, a great deal of 
^ ■ effort has been devoted to implementing Weinberg's original prescription [l| for the problems of few- 

nucleon systems, with the two-nucleon system as the starting point j^-Q] (for more general reviews, 
see Refs. [9- 14]). While the power counting of pion exchange diagrams follows the paradigm of 
chiral perturbation theory (ChPT), i.e., chiral EFT in the single- nucleon sector, estimating the 
size of NN contact interactions often requires assumptions beyond chiral symmetry. 

Assumed in Weinberg's power counting (WPC) is what we refer to as naive dimensional analysis 
(NDA): each derivative on or each power of pion mass dependence of the Lagrangian terms is always 
suppressed by the underlying scale of chiral EFT, My ~ m a , where m q is the mass of the a meson. 



While plausible, this assumption was questioned in a number of works 



1514231] , and was shown to be 
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only partially correct from the perspective of renormalization group (RG) invariance; WPC does 
not have enough NN contact interactions, or counterterms, to renormalize the nucleon-nucleon 
(NN) scattering amplitudes at a given order. One mechanism to spoil the renormalizability of 
WPC is the singular — diverging at least as fast as 1/r 2 — attraction of the tensor force of one-pion 
exchange (OPE): — 1/r 3 at r — > 0. In S and P waves, the triplet channels subject to this singular 



23|, 



attraction include one uncoupled, 3 P , and two coupled, 3 S, - 3 D, and 3 P 2 - 3 F 2 . In Ref. 
we used 3 Pq to investigate the modification to WPC under the guidance of RG invariance. In 
this paper, we report a complete study of the triplet channels, in which we continue our efforts 
to modify WPC at the subleading orders, in a generalization of Ref. [24]. Interestingly, we reach 



conclusions that differ in some aspects from a parallel investi gati on in Refs 



■3,i 



2J. 



25-231, it is well accepted that in S 



1<|. Even 



Except for the attempts to treat OPE as perturbation 
and P waves the leading order (LO) amplitude requires the full iteration of OPE [l£ 
though the nonperturbative unitarity requires any nonperturbative, nonrelativistic T-matrix to 
scale as Q , we choose to label the LO as 0(1) so that one does not need to change the standard 
ChPT notation for power counting pion exchange diagrams. What is more consequential is that 
we denote different orders of the EFT expansion by its relative correction to the LO, i.e., the 
next-todeading order (NLO) by 0{Q /M^) or 0(Q) for short, and next-to-next-todeading order 
(NNLO) by 0(Q 2 or 0(Q 2 ), and so on. 1 RG invariance, or, more specifically, invariance of 
the amplitude with respect to the ultraviolet (UV) momentum cutoff inherent in the Lippmann- 
Schwinger equation, demands as well one counterterm to be fully iterated in the partial waves 
where the tensor OPE is attractive 18j |. In P waves, these singular attractive channels include 
3 -Po an d 3 P 2 — 3 -^2- However, WPC considers the counterterms in these channels to be subleading, 
for they are second-order polynomials in momenta. Stated differently, the leading counterterms in 
3 P and 3 P 2 - 3 F 2 are underestimated in WPC, and RG invariance requires them to be enhanced 
by 0(M^j/M|q). Here, M\ a refers to a cluster of infrared mass scales that include the pion decay 
constant f n ~ 92 MeV, the pion mass m n ~ 140 MeV, and certain combinations of them. 

Although not one of the centerpieces of WPC, the indiscriminate, full iteration of different 
order potentials as a whole has been the standard practice in its implementations. However, the 
ordering of the potentials according to their matrix elements for small momenta is not always 
valid in the nonperturbative treatment in which the intermediate states could reach My, where 



1 Here NLO and NNLO are defined differently from a more conventional notation 0, where NLO is 0(Q 2 ) and 
NNLO is 0(Q 3 ). 
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the higher-order potentials usually have larger matrix elements. In the nonperturbative setup, it 
seems to take intricate cancellations for higher-order potentials to eventually contribute less to 
the low-energy, on-shell amplitude. But it is far beyond the scope of our paper to decide whether 
or under what conditions these cancellations will happen. We refer readers to Refs. 



28, 



2g for 



discussion regarding renormalization and power counting in the nonperturbative treatment. 

To minimize the interference between lower- and higher-order potentials in the UV region, 
we choose the natural way to go beyond the LO, that is, to treat the sub leading interactions as 
perturbations. Now that the potentials from different orders are no longer on an equal footing, it is, 
as we will see, much easier to separate in the UV region the contributions of higher-order potentials 
from those of lower-order ones. If the subleading interactions are too strong to be perturbative, they 
are simply not subleading in a bona fide EFT. Reference |2J] explained the perturbative formalism 



with a toy model: — 1/r 2 as LO and ±l/r 4 as 0(Q 2 ) long-range potentials. More importantly, the 
general lesson drawn from the study of Ref. [2jy] , referred to in the paper as modified NDA ( NDA ) , 
is that in the case of the LO long-range potential being singular and attractive, the subleading 
counterterms (SCTs) are enhanced relative to NDA by the same amount as the LO counterterms; 
as the long-range force gets an 0(Q 2 ) correction, so do the contact operators that have two more 
derivatives than the LO counterpart. 



The validity of NDA is confirmed by renormalization of uncoupled 3 Po up to 0(Q 3 ) |22j, l23l ]. 
which shows that the leading long-range potential being —1/r 2 is not essential to the applicability 
of NDA. Now 3 Si — 3 D\ and 3 i-2 — 3 -^2 pose interesting questions as to the extension of NDA to 
the coupled-channel problems. Take 3 Si — 3 D\ as an example. The LO counterterm is a constant, 



C-is 1 , and, according to NDA, two second-derivative terms will turn up at 0(Q 2 ). Therefore, NDA 
suggests a total of three counterterms up to 0{Q 3 ) in each of 3 Si — 3 D± and 3 i-2 — 3 ^2 ? which will 
be discussed in more details in Sec. IIII CI However, using a coordinate space setup, Refs. 21, zj] 
concluded that there must be six counterterms in 3 Si — 3 D\ or 3 P2 — 3 -^2 up to 0(Q 3 ). If this 
proliferation of counterterms in the coupled channels is true, then the predictive power of nuclear 
EFT is further weakened. 

We carry out in the paper a momentum-space calculation to verify the power counting based 



on NDA. In particular, we are interested to see which can be confirmed in the coupled channels: 



the proliferation of six counterterms or three counterterms inferred from NDA. Our principle of 
establishing power counting is summarized as follows: 

(i) The size of pion exchanges is decided by the non-analytic part of the corresponding Feynman 
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diagram, which is rightly captured by WPC. 

(ii) We promote counterterms over WPC only if RG invariance requires it. That is, when a 
counterterm is not needed for renormalization, its counting will follow NDA. 

The rationale for the second point is that RG analysis in terms of the floating momentum cutoff 
of the NN intermediate states touches upon only (nonrelativistic) nucleon momenta and it does 
not "know" anything about the contributions of heavy mesons that are integrated out in the first 
place [30! ] . 

A study of power counting of chiral AW forces is not complete without the singlet channels. 
But the drastically different short-range behavior of OPE in the singlet (1/r) and triplet (1/r 3 ) 



channels signa 
further study 



s different structures of counterterms. Therefore, we leave the singlet channels to a 
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We also leave out D and higher waves in our analysis, except for 3 Lq and 3 F 2 which are, 
respectively, coupled to 3 S± and 3 P2- The reasons are as follows. First, the conceptual issues 
concerning renormalization and power counting can be well illustrated by S and P waves. Second, 
it is debatable whether or to what extent OPE is perturbative in D waves 
answering this question is beyond the scope of our paper. 

After briefly reviewing the LO and establishing our notation in Sec. UH we will establish in 
Sec. imthe SCTs in 3 Pi, 3 P 0) 3 Si — 3 Di and 3 P2 — 3 -^2 by examining the cutoff dependence of the 
subleading amplitudes. Finally, we offer a discussion and a conclusion in Sec. IIVI 

II. LEADING ORDER 

OPE is the leading long-range NN interaction, 

V L ( q ) = V 1 M = -j f ,r 1 .T 2 _, 2 + r?4 , (1) 

where q = p' — p is the difference between the outgoing [p') and the incoming (p) momenta in 
the center-of-mass frame, the axial vector coupling constant qa = 1.29, the pion decay constant 
f n = 92.4 MeV, and the pion mass = 138 MeV. Its coordinate space version will be useful in 
the discussion, 

Vi^r) = Kn • t 2 [T(r)Si2 + Y(r)a x ■ a 2 ] , (2) 



5 



where 



™% f 9a 



127T \,4/2 



A ^ - 777 ( 772 ) ' ( 3 ) 



T(r) 



„—m n r 



m^r (m v r) 2 



(4) 



y( r ) = , (5) 

and 

5i2 = 3(<ri • r)((T2 • r) - <Ti • <t 2 . (6) 

The tensor force T(r) has an inverse cubic short-range core, 1/r 3 , but it contributes to only triplet 
channels. 

In the lower partial waves, the LO amplitude is obtained by the full iteration of OPE 
and necessary counterterms, through solving the Lippmann-Schwinger equation. For the coupled 
channel with total angular momentum j, the off-shell Lippmann-Schwinger equation reads 

2#V,p; *) = VS\ P ',P) + 2 ~m N £ / A «k . 2 -) f l } K ff ie > ( 7 ) 

with Z, Z', and Z" running over j — 1 and j + 1, fe as the center-of-mass momentum, and A as 
the momentum cutoff. Extension to the uncoupled channels is straightforward. The product of 
the Schrodinger propagator and the integral measure scales as m^Q. Since OPE scales, more or 
less casually, as (mTv^io) -1 , with M\ Q a certain combination of and mj\r, OPE must become 
nonperturbative when Q ~ M\ Q . 

Depending on the sign of the matrix element of S12, the OPE tensor force drives NN contact 
interactions in very different ways 

{lsj\Si2\lsj) is positive, where Z is the orbital angular momentum, s = 1 is the spin, and j is the total 
angular momentum, the OPE tensor force is ~ +l/r 3 . If one picks up the regular solution to this 
repulsive potential, the wave function dies off exponentially near the origin. As a consequence, the 
sensitivity to the UV cutoff vanishes very quickly; thus, there is no need for an extra counterterm 
to absorb the cutoff dependence jq|. This is in agreement with WPC because, according to WPC, 
the first counterterm in 3 Pi — the lowest repulsive, uncoupled triplet channel — appears at 0(Q 2 ). 

When (Isj \ S12 \ Isj) is negative, the OPE tensor force overpowers the kinetic energy and the 
centrifugal barrier, causing the NN system to "collapse" The mathematical origin of this 

pathology is the simultaneous existence of two equally good solutions to the Schrodinger equa- 
tion [a], which eventually lead to ambiguity in predicting physical observables. Or, in terms of the 
Lippmann-Schwinger equation, the scattering amplitude is very sensitive to the UV cutoff. 



The modern-day interpretation of renormalization offers a cure to this sort of pathological 
potential with singular attraction near the origin: supplementing short-range interactions Vs rather 
than naively extrapolating the long-range interaction to short distances. The model-independent 
treatment involves arranging Vs as counterterms that run with the UV cutoff A in such a way that 
physical observables do not depend on A 



16, 



19, 



35 



36| . Although this is not a complete innovation 



in the context of chiral EFT since NN contact terms are always part of the chiral Lagrangian, the 
consideration of RG invariance offers a priori insights into gauging the importance of counterterms 
at a given order. 

In 3 Po where I > 0, the singular attraction of OPE requires at 0(1) a counterterm, which is, 
however, considered 0(Q 2 ) in WPC for it is a second order polynomial in momenta [181 ]. 

( 3 Po\v^\ 3 p ) = a PoP 'p, (8) 

where p' (p) is the magnitude of p' (p). Now that Czp^p'p is at LO with OPE ~ 47r/(mjv-Mio)) 
C3p must scale as 

where Ait/tun is introduced to cancel a common factor that usually concurs with loop integrals 
involving NN intermediate states. For comparison, WPC considers C.ip p'p to be of the same 
order as the leading two-pion exchange, which in turn is counted as an 

C(Q 2 / M hi) correction to 

OPE; therefore, WPC prescribes Q] 

In the coupled channels, it is convenient to write the tensor OPE in coordinate space as a 2 x 2 
matrix in the basis of two coupled orbital angular momentum states, I = j ± 1 js?]]. In 3 Si — 3 Z?i 



and 3 P 2 - 3 F 2 , 



-6^2 , 

V T in( 3 Si - 3 Z>i) = K | T(r), (11) 

-6y/2 6 



W 3 P 2 - 3 F 2 ) = ^H 2 " \ T (r). (12) 



5 



-8 6^ 



There is a domain of r near the origin in which T(r) dominates over the centrifugal barrier and in 
which diagonalizing Vtx-k also diagonalizes the Schrodinger equation. OPE matrix elements in all 
of the coupled channels share one property: there always is one attractive and one repulsive eigen 
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subchannel. Therefore, one needs to summon at LO a short-range input to counter the singular 
attraction in one of the subchannels. This accords with WPC in 3 Si — 3 D\ but calls for amendment 
in 3 P 2 - 3 F 2 Q- 

Although the counterterms are easily formulated as polynomials in momentum space, the above- 
mentioned diagonalization cannot be trivially realized therein. The counterterms in 3 Si — 3 D\ and 
3 P 2 — 3 F 2 have the generic momentum space form 

, q n q / CW + £?3 Sl (p' 2 + p 2 ) EsdP 2 \ 

( 3 S 1 - 3 D 1 \V S \ 3 S 1 - 3 D 1 )=[ Sl bD/ l +.., (13) 

\ E SD p' 2 Fs Dl p' 2 p 2 J 

, l Cz P „ + Dz Po ip' 2 + p 2 ) E PF p 2 
( 3 P2- i F 2 \V s \ 3 P 2 - 3 F 2 )=p'p[ Pa ^ ^ PF ^J+.... 

s f 2 : 



£ PF p' 2 F 3F2 p' 2 p 2 



As shown in Ref. it is not necessary to design a counterterm that exclusively acts on the 
attractive subchannel of Vrin dUJi the C term alone will properly renormalize all of the T-matrix 
elements. 

While the C term in 3 Si — 3 D± does not violate WPC, the promotion of the C term in 3 P 2 — 3 F 2 

( 3 S 1 - 3 D 1 \V^\ 3 S 1 - 3 D 1 )=[ C ^ °| , (15) 

0, 



asks for an enhancement of 0{M?JM?^: 



r (0) n - 

( 3 P2- 3 F 2 \V^\ 3 P 2 - 3 F 2 )=p'p[ 3P2 | . (16) 

0, 



Stated differently, the renormalized Cz Sl and Cz p 2 scale as 



C5S ^^Wo> aP2 ~^4' (17) 



whereas WPC differs for Czp 2 , 



^WPC 4vr 1 

^ ~^^r (18) 



III. SUBLEADING ORDERS 
A. Generalities 



It has been long known that the two-pion exchanges (TPEs) with chiral index v = vertices 
(TPEO) give an 0(Q 2 ) correction to OPE, and that the TPEs with one insertion of v = 1 vertices 
(TPE1) lead to 0(Q 3 ) long-range potentials. TPEs are not uniquely defined because different 
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381 ] — applied 



regularization schemes — such as dimensional and spectral function regularizations 
to two-pion-exchange Feynman diagrams may lead to different expressions. To remove the arbi- 
trariness in the choice of regulator, one should always pair TPE expressions with a momentum 
polynomial. We refer to this polynomial as the "primordial" counterterm for TPEs. By definition, 
the primordial counterterm is power counted as the same order as the corresponding TPEs, which 
is exactly the content of WPC. 

In addition to the general rationale given at the end of Sec. HI the concept of primordial coun- 
terterms serves as a complementary argument against demoting counterterms in the cases where 
they are not needed for renormalization in the context of the Lippmann-Schwinger equation. 

TPEs, as the two-pion-exchange diagrams evaluated in the plane-wave basis, are not the only 
sources driving contact interactions when high momentum modes are integrated out. There are 
two other classes of diagrams that contribute to the evolution of the SCTs. One is insertions of 
TPEs into the LO T-matrix, which, after being properly renormalized, generates an 0(Q 2 ) or 
0(Q^) correction to the LO [see Eq. (|19p ]. The structure of divergence in these diagrams has been 



used as the primary tool to gauge the SCTs 



21 



24|. 



The other mechanism driving the SCTs is the LO amplitude itself. With the necessary coun- 
terterms, if any, the cutoff dependence of the LO amplitude vanishes at A — > oo; there normally is 
residual cutoff dependence at finite As. When A is rescaled to a smaller value, A', the integrated- 
out momentum modes in the interval of (A', A) will, for the most part, contribute to the evolution 
of the leading counterterm, but a perfect RG invariance would request the SCTs to evolve as well, 
so that there would not even be a small cutoff dependence. References Jj], [^(J studied this mech- 
anism using Wilson's RG equation, but their assumption about the fixed-point solutions to the 
RG equation seems to be at odds with the limit-cycle-like behavior of the LO counterterms of the 
attractive channels 

A simpler way to gauge the LO-amplitude-induced SCTs is to consider how the LO residual 
cutoff dependence scales against A. According to Ref. the LO residual cutoff dependence 
in the attractive triplet channels is 0(A~ 5 / 2 ), which means our ignorance of the LO-amplitude- 
induced SCTs is smaller than 0(Q 2 /M?A — the corrections brought by TPEO and its primordial 
counterterms. As a consequence, there is no need to have a nonvanishing 0{Q) counterterm that 
is not accompanied by any long-range force. 

Yet another evidence of vanishing 0(Q), though a posteriori, is that the LO error of the 
3 Si — s Di mixing angle scales as k 2 instead of k, \ef FT — e^ WA | oc k 2 , as seen in Fig. 6 of Ref. \v\ . 

In conclusion, the counting of the SCTs in the triplet channels will be decided by the larger one 



9 



of (i) their primordial size 0(Q 2 ), and (ii) what the divergence of one insertion of TPE requires. 

With 0{Q) vanishing, the on-shell <D(Q 2 ) and 0(Q 3 ) T-matrices, and , are calculated, 
respectively, by one insertion of 0(Q 2 ) and 0(Q 3 ) potentials into , 

T^' 3 ) (jfc, jfe) = V < 2 > 3 ) (k, k) + -m N [ A dn k 2 3 ) (jfe, k) J^^'^ 

7T J k z -K z + ie , , 

+ —m 2 N / dKdit'K 2 *,' 72 l^-y^' 3 ) — v ' y , , 

^ J J hr - k z +ie k 2 - k' z + ^e 

where T-matrices and Vs are understood as 2 x 2 matrices for the coupled channels. This is of 
course nothing more than the first-order distorted-wave expansion. 

Treating the subleading potentials as perturbations will no doubt break the exact unitarity 
of the 5-matrix, as does any perturbation-theory-based calculation. But in a consistent power 
counting scheme, the violation of unitarity is of higher order. This makes it slightly nontrivial to 
extract the phase shifts and the mixing angles from the expanded EFT T-matrix. Although it has 
been covered in the literature, we list in the Appendix [A] the useful formulas for convenience of 
reference. 



B. Uncoupled Channels: 3 Pi and 3 Pq 



In 3 Pi, OPE (with short-distance behavior +l/r 3 ) leads to an LO wave function exponentially 
suppressed near the origin: ~ exp[— (or) 1 / 2 ], with a as a positive mass scale [^(J. The exponential 
damping of the LO wave function would eliminate the singularity of TPEO (~ 1/r 5 ) and TPE1 
(~ 1/V 6 ), even without any counterterm. But, as we argued, we will not demote any counterterm 
with respect to WPC. Therefore, the 0(Q 2 ) and 0{Q i ) 3 Pi counterterms are 



< 3 ^i|V". 



(2, 3) | 3 
S 



a 



(0, 1) 



3 Pl VP- 



(20) 



The splitting of C , 3p 1 into different orders does not mean that we will take more than one input for 
C f 3p 1 ; it only reflects the possibility that the "bare" values of the counterterms could be modified 
by the short-range core of TPEs. 

Perturbative renormalization at subleading orders in 3 Po was first studied in Ref. 
a parallel work of ours 



channels as Ref. 



221 ] and in 



231 ] . Although we have reached the same conclusion about the uncoupled 
22|, we include here, for completeness, our analysis of 3 Pq 231 ] . 



After the LO amplitude is renormalized with C3p , the LO wave function in 3 Pq can be approx 



imated near the origin in powers of k 2 up to a normalization factor 



13, 



39, 



4l( 



V>f (r) 



AV 

r 



u + k 2 r 2 Jj Ul + O(k A ) 



(21) 
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where A = ^^^jz > u o an d u\ are oscillatory functions in terms of r / A and <p with amplitudes ~ 1 , 
and 4> is the phase between the two independent solutions and is related to Cs,p . Combined with 
the short-range behavior of TPEO, ~ 1/r 5 , and TPE1, ~ 1/r 6 , the superficial divergence 



of one insertion of TPE is estimated on a dimensional ground 



22. 



23j: 



<3P = ^ (0) l4V 0) >3p ~ / drr 2 \^(r)\ 2 \ 

~ ao(A)A 5/2 + f3o(A)k 2 + 0(fc 4 A" 5 / 2 ), (22) 
^ = ^ (0) I^V {0) >3P ~ l iA drr^\r)\>± 

~ «i(A)A 7 / 2 + /3i(A)AA: 2 + 0(A; 4 A- 3 / 2 ) , (23) 

where Qo,i(A) and /?o,i(A) are oscillatory functions diverging slower than A. 

The presence of two divergent terms suggests that (i) the running of the LO counterterm Czp Q (A) 
needs to be corrected at higher orders, and (ii) one SCT, D<,p p'p{p' 2 + p 2 ), needs to be enlisted. 
The fact that Dz P arises at the same order as TPEO leads us to arrange counterterms at 0(Q 2 ) 
and 0(Q 3 ) as follows: 

( 3 Po|# 3) I 3 Po) = C? p f p'p + p'p(p' 2 + V 2 ) ■ (24) 

We see that the enhancement of Dz Pq is the same as that of C3p Q : 0(M 2 Yi /M 2 Q ). 



The lesson of power counting learned here coincides with the conclusion of Ref. 



2J]. Although 



NDA fails to prescribe a counterterm at LO, we could use NDA to determine how the SCTs scale 
when subleading long-range potentials are taken into account, which states that the enhancement 
of each SCT is the same as the LO counterterm. 



C. Coupled channels: A S-y i D l and A P 2 A F 2 

In the coupled channels, the LO wave functions are dominated at short distances by the at- 
tractive subchannel. Because there are three independent on-shell T-matrix elements, calculating 
the superficial divergence of TPEs using Eqs. (|22|) and (|23|) gi ves rise to six divergent terms, with 
two for each T-matrix element, as shown in detail in Ref. 22 1. Reference [3] proposes to use six 



counterterms to cancel these divergent pieces on a one-to-one basis. Although it guarantees RG 



invariance, lost is the regularity of power counting enjoyed by NDA in the uncoupled channels. 
On the other hand, the presence of six divergent terms does not necessarily mean that one 



must have six counterterms to achieve RG invariance. NDA suggests that when the LO long-range 
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potential gets 0(Q 2 ) correction going from OPE to TPEO, so do the SCTs; therefore, the 0(Q 2 ) 
SCTs in the coupled channels should be the D and E terms in Eqs. (I13p and (I14p . which have two 
more derivatives than the LO counterterm. We propose the following power counting based on 



NDA for the SCTs of the coupled channels: 
(i) In 3 Si — 3 Di, we do not change WPC, 



( 3 5 X - ^v^sx - 3 Di> = I C ^ 3) + D ;° s \ ){p ' 2 +p2) E ™ )p2 

4°dV J (25) 



(ii) In 3 P2 — 3 F2, an enhancement of 0{M 2 i /M 2 ) leads to 
( 3 P 2 - 3 F 2 \V^ 3) \ 3 P 2 - 3 F 2 )=p'p 



3d 3p|T/(2,3)|3p 3 pU „'n I 3 ^ + 3 P 2 ^ + P ) E pF p 



4fV / (26) 



An analytical proof of renormalizability with the above counterterms is difficult because the closed 
form of the LO T-matrix is not available, so we will resort to numerical experiments in Sec. IIII Dl 
to test RG invariance, or the lack thereof. 

At a given energy, there are three scattering parameters — two phase shifts and one mixing 
angle — to be extracted from the 2x2 T-matrix (|Alj) . With determined by three inputs from 
partial- wave analysis (PWA) and Eq. (|A6]), Eq. (|19p provides a group of linear equations to solve 
for C( 2 \ T)(°), and E^°\ Though not obvious, it is straightforward to show that the three linear 
equations built from the three PWA inputs at the same energy are linearly dependent. Therefore, 
the needed three PWA inputs must be incorporated from at least two different energies. Not 
surprisingly, the same also applies to T^ 3 \ 



D. Numerics 



42| in the literature with slight differences in how double 



We will use TPEs without the explicit delta-isobar to demonstrate renormalization. There are 
a few versions of TPEs 

counting is avoided 43 |. For definitiveness, we use the version in Ref. jjj, i.e., delta-less TPE 
expressions with dimensional regularization. We adopt the following low-ener gy constants for the 
v = 1 ttttNN seag ull couplings (GeV- 1 ): c\ = -0.81, c 3 = -4.7, and c 4 = 3.4 [4J. With the delta 
integrated out, we expect the EFT expansion to break down around Q ~ 5, where 5 ~ 300 MeV 
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FIG. 1: (Color online) With the SCTs (J3DJ), the 3 P 1 phase shifts (a) as a function of T lab up to 0(Q 3 ), and 
(b) as a function of A for Ti a b = 100 MeV. In the legend of (a), the number following the symbols is the 
cutoff value in GeV. The red dots are from the Nijmegen PWA. 



is the delta- nucleoli mass splitting. Such a breakdown scale is noticed in, e.g., Ref. [45J through 
"deconstruction." 

Figure Q] shows the 3 P\ phase shifts as a function of laboratory kinetic energy Tj a b for two 
different cutoffs and a function of the cutoff at Ti a b = 100 MeV. The sharp momentum cutoff is 
chosen as the regulator throughout the paper: 6 (A — k), where K is the magnitude of the loop 
momentum, as defined in Eq. ([7]). The values of Cgp^ are solved for by a fit of the phase shift 



to the Nijmegen PWA 



461 ] for Ti a b = 50 MeV. The cutoff independence is not surprising since the 
LO wave function is exponentially suppressed at short distances. The large shift from 0{Q 2 ) to 
0(Q S ) indicates the influence of the uncertainties of itttNN coupling constants q. Nevertheless, a 
decent agreement with the PWA up to = 100 MeV is obtained. 

3 Pq has been thoroughly studied in our momentum-space framework in Ref. [23], which con- 



firmed the RG invariance of NDA for the uncoupled channels; see (|24p . We refer the reader to 



Ref. 



231 ] for the numerical results. 



Now we move on to the coupled channels. Figure [2] shows the phase shifts of 3 <Si — 3 -Di and 
the mixing angle e\ as functions of at 0(Q 2 ) and 0(Q 3 ). The values of the counterterms are 
determined such that the EFT curves reproduce the Nijmegen PWA for 63$ 1 at Tj a b = 30 and 
50 MeV and e\ at 50 MeV. Consequently, the 3 D± EFT phase shifts are predictions. The first 
indication of the cutoff independence is the closeness of two EFT curves with A = 1.5 and 2.5 
GeV. The plot of ei shows a larger cutoff dependence toward higher energies, but the fact that 
the A = 2.0 GeV curve is closer to A = 2.5 GeV than A = 1.5 GeV suggests that the cutoff 
independence is finally achieved at larger As. 
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FIG. 2: (Color online) With the SCTs (f25j) , the 3 5i, 3 £>i phase shifts and the mixing angle t\ as functions 
of Ti a b at 0(Q 2 ) and 0(Q 3 ). In the legend, the number in front of the line symbols is the cutoff value in 
GcV. 

The cutoff independence is more clearly demonstrated in Fig. [3l which shows the phase shifts 
and the mixing angle at Ti a b = 40 and 100 MeV as functions of A. The residual cutoff dependence 
is still visible at lower As, but it is much smaller than the size of the corresponding EFT correction. 
That is, the corrections are meaningful even at the lower cutoffs because they are not washed out 
by the cutoff uncertainties. 

Figures H] and [5] show the 3 P2 — 3 -?2 phase shifts and the mixing angle e 2 at 0(Q 2 ) and C(Q 3 ) as 
functions of T\ a b and the cutoff, respectively, with the SCTs (|26p . Similar to the case of 3 S*i — 3 -Di, 
we fit (5 3 p 2 at T lab = 30 and 50 MeV and e 2 at 50 MeV to the PWA values. 

The main goal of this paper is to verify the RG invariance of our power counting in the UV region, 
so the cutoff window was chosen such that the EFT curves start to show the cutoff independence. 
For smaller cutoffs not shown in the plots (1.2 > A > 0.6 GeV), the general trend is similar to 3 Po 
(Fig. 2 of Ref. 23J): the 0(Q 2 ) EFT curve is the first to become cutoff independent, while the LO 
is the latest. 

A good fit to the PWA up to 100 MeV is presented at 0(Q 3 ) in 3 Si - 3 D X . Without special 
effort to improve the fits, the EFT result agrees less well with the PWA in 3 P 2 — 3 ^2- We think 
that this is largely owing to a disappointing LO, which departs quickly from the PWA as the energy 
increases. The unusually small 3 P 2 scattering volume, azp 2 ~ —0.28 fm 3 , compared with a of other 
P waves, \a\ ~ 1.5-2.8 fm 3 H, 

is suggestive of a certain amount of fine tuning, which calls for 
a more sophisticated fitting strategy. 

To assess the feasibility of improving the fitting quality, we refit the counterterms to the PWA 
inputs at higher energies: <5 3 p 2 at 50 and 100 MeV, and e 2 at 50 MeV. The updated 3 P 2 - 3 F 2 EFT 
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FIG. 3: (Color online) With the SCTs ([25]), the 3 Si, 3 Di phase shifts and the mixing angle ei at T lah = 40 
and 100 MeV, as functions of the momentum cutoff. The dashed, dot-dashed and solid lines are 0(1), 0(Q 2 ) 
and 0(Q 3 ), respectively. 
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FIG. 4: (Color online) With the SCTs (|26|) . the 3 P2, 3 F 2 phase shifts and the mixing angle e 2 as functions 
of Tiab at 0(Q 2 ) and 0(Q 3 ). The symbols are explained in the legend of Fig. [5] 
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FIG. 5: (Color online) With SCTs (|26|). the 3 P 2 , 3 F 2 phase shifts and the mixing angle e 2 at Ti ab = 40 and 
100 MeV, as functions of the momentum cutoff. The symbols are explained in the caption of Fig. [3] 

phases are plotted in Fig. [6l Although the EFT convergence still breaks down at lower energies 
than in 3 Si - 3 L>i, a good fit to the PWA until 130 MeV is achieved at 0(Q 2 ) and C(Q 3 ), and it 
is comparable to the WPC-based calculation with the same TPEs, which is shown in Ref. Q|. We 
notice another fitting strategy used in Ref. 



221 ]. which sacrifices the LO near threshold in order 
to facilitate better agreements with the PWA at higher orders. This amounts to tuning the LO 
counterterm to further reduce the attraction of OPE. 

Overall, the breakdown scale implied in the numerical results is consistent with our expectation 
for the delta- less theory, k ~ 5 ~ 300 MeV, with the exceptions of 3 Pi and 3 p2 — 3 ^2- Even for 
these two channels, one cannot help wondering whether the delta can bring some attraction from 
0{Q?) to 0{Q 2 ) and improve the convergence of EFT expansion 48 1. 
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FIG. 6: (Color online) With the SCTs (|26[) . the 3 P2, 3 -p2 phase shifts and the mixing angle e 2 as functions 
of Xiab with A = 1.5 GeV. The values of the counterterms are adjusted to reproduce the PWA values for 
5a p 2 at Ti a b = 50, and 100 MeV and e 2 at 50 MeV. The dots arc from the Nijmcgcn PWA. The dashed, 
dot-dashed and solid lines are 0(1), 0(Q 2 ), and 0(Q 3 ), respectively. 

IV. DISCUSSION AND CONCLUSION 



We have studied up to 0(Q 3 ) the structure of counterterms of chiral NN contact interactions in 
the triplet channels, with S and P waves as the examples. The essential guideline we have followed 
is to promote counterterms over WPC when RG invariance requires it. 

We found that the scaling of SCTs are mainly driven by the interplay between TPEs and the 
LO wave functions. A direct consequence is that 0(Q) by contact interactions alone vanishes. The 
resulting arrangement of the counterterms in the studied channels are given by Eqs. (|20h . (I24p . 



(|25p . and (|26p . which can be very nicely summarized by NDA: the SCTs are enhanced by the same 
amount as the LO counterterm so that the whole tower of counterterms with the same quantum 
number is shifted uniformly. While this means that WPC remains intact in 3 Si — 3 -Di and 3 Pi, it 
requires an enhancement of 0(M&/M?) to all counterterms in 3 Po an d 3 i"2 — 3 -p2- 

It is interesting to compare the chiral AW forces with the three-body system that has only 
contact interactions. When the two-body 5"-wave scattering length a<i — > 00, the three-body "pion- 



less" theory can be mapped onto a dual two-body theory with — 1/r 2 lon g-ra nge force anc 



interactions that represent three-body operators in the original system 



49( . Reference 



contact 



pro- 



posed a power counting similar to NDA for the three-body contact interactions. However, the 



resemblance between chiral EFT forces and the system investigated in Ref. [50(] is not perfect be- 
cause the long-range interactions beyond the leading —1/r 2 in the dual two-body system, if any, 
are resummed nonperturbatively instead of being treated as perturbations. Therefore, it is not 
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50l | is driven by the LO or subleading 



clear to us whether the scaling of SCTs obtained in Ref. 
long-range interactions. 

Our finding that three counterterms are needed in the coupled channels up to 0{Q i ) differs 



from that of Refs. 



21 



, |22| , which also adopted the perturbative approach on top of the nonpertur- 



bative LO, but concluded instead that six counterterms are necessary for renormalization purpose. 



(Although it is speculated in Ref. 22| that it might be possible to reduce the number of short- 
range parameters, an alternative is not offered unless the higher-wave component is treated in 
perturbation theory.) However, that there are only two second-derivative terms — Dzg and E$n 
in (|T3|) — suggests that one could correlate these six divergent pieces in a model-independent way, 
which is justified by the numerical evidence of RG invariance shown in Figs. [3]and[5l 

Without a dedicated effort to fine tune the fits, our results show a good agreement with the 
Nijmegen PWA up to Ti a b ~ 100 MeV. Regardless of the comparison with the PWA, the relatively 
large deviation from 0(Q 2 ) to 0(Q 3 ) in 3 Pi and 3 i-2 encourages one to hope that a delta-ful EFT 
can improve the convergence by including the delta-isobar as explicit degrees of freedom. Aside 
from the debatable issues of power counting counterterms, the delta-ful nuclear forces have been 

, 5l| and, on a more qualitative 



shown to achieve a more rapid convergence in the two-nucleon 
level, the three-nucleon 



29 



521 ] sectors. 

Although there have been many efforts to derive the delta-ful TPEs 
desirable to update the extraction of low-energy constants in the delta-fu 



2- 



33 



531 ]. it appears 



chiral Lagrangian from 



54. 



5511 where the effects 



nN scattering through the chiral EFT description around the delta peak 
of the delta are most prominent. In addition, the formulation of the delta-ful chiral Lagrangian 
may need to be reexamined in light of the discussion in Ref. [561 ]. 
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Appendix A: Perturbative Relation between Phase Shifts and S'-matrix 



In distorted- wave expansion, the unitarity of the S'-matrix no longer rigorously holds. Therefore, 
the Stapp parametrization (571 ] of the S'-matrix in the coupled channel, 

cos(2e)e 2i(Sl i wx(2e)e^ 5l+s ^ 
ism(2e)e^ 5l+s ^ cos(2e)e 2i<52 



S 



(Al) 



needs to be adjusted, although it is still valid at 0(1) where the LO potential is fully iterated. 
Here, 5\ and 82 are the phase shifts of the partial wave with I = j — 1 and I = j + 1, respectively, and 
e is the mixing angle. Suppose that with S*- 1 ^ vanishing, the S-matrix and the phase parameters 
have the following expansion (labeled by the usual EFT order): 



S = S (0) + 5 (2) + 5 (3) + 



'1,2 



,(o) ■ A 2 ) 



+ + 



(A2) 
(A3) 
(A4) 



Expanding both sides of Eq. (|A1[) . one finds that 



o(2) o(2) N 
D ll D 12 



c(2) o(2) 
°21 °22 



(A5) 



with 



0(2) 
^11 



J2i5 



(0) 



5 



(2) 
12 



? (2) 
'22 



(0) 



-2e< 2 > sin(2e( ))+i2,5[ 2) cos(2 e (°); 

2e( 2 ) cos(2e(°)) + i(5® + jf) sin(2eW)l , (A6) 
-2e( 2 )sin(2e(°))+i24 2) cos(2eW)" . 

Replacing the superscript ^ with ( 3 ) in Eqs. (lA5j) and (|A6p . one obtains the relations for 0(Q 3 ). 
To convert the T-matrix to the S-matrix, notice the normalization adopted in the paper: 



S = 1 - 2ikm N (tW+T® + 



(A7) 
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